Group theoretical analysis of double acceptors in a magnetic field: identification of 

the Si:B + ground state 
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A boron impurity in silicon binding an extra hole is known to have only one bound state at an 
energy of just below 2 meV. The nature of the Si:B + ground state is however not well established. 
We qualitatively analyze the behavior in a magnetic field of isolated acceptors in a tetrahedral lattice 
binding two holes using group theory. Applying these results, we analyze recent measurements and 
conclude that the ground state of B + is most compatible with a non-degenerate Ti state. 

PACS numbers: 71.70.Ej,03.65.Fd,71.55.Cn,71.70.Ch 



INTRODUCTION 

A neutral boron acceptor in silicon is able to weakly 
bind an extra hole, resulting in a positively charged ion 
(B + ). This entity is an example of a positively charged 
acceptor, commonly denoted by A + , which is the counter 
part of the better known negatively charged donor D - . 
Both are related to their atomic-physics analogue, the 
negative hydrogen ion H~ . The energy states associ- 
ated with these ions are very shallow and spatially large. 
When their concentration is sufficiently high, their over- 
lapping wave functions can form an upper Hubbard band 
flj and they play an important role in electronic trans- 
port in semiconductors at low temperatures. Since a few 
years, electronic states of individual dopant atoms gained 
renewed interest due to their prospective application in 
Si-based solid state quantum computing 2] . 

Neither theoretically nor experimentally much work 
has been done on the B + -state. Optical spectroscopy 
is difficult due to the small ionization energy (less than 
2 meV 3]). In particular, the nature of its ground state 
is not well-established and to our knowledge no results 
have been published on the magnetic field dependence of 
the energy levels. 

The purpose of this paper is twofold. First, we present 
a general group-theoretical study of the magnetic field 
dependence of two-hole states in tetrahedral semiconduc- 
tors. To our knowledge, such an analysis has not been 
published before. Second, because our analysis includes 
all possibilities for the B + ground state, it enables us to 
compare our results with our recently published measure- 
ments of B + in a magnetic field Q and to draw conclu- 
sions about the nature of the B + ground state. 



BACKGROUND 

The nature of the energy levels of a neutral boron ac- 
ceptor (B°) in silicon is well-known and the Zeeman 
effect in B° has been studied in detail, both theoretically 
and experimentally Q. The B-impurity is located at 
substitutional sites of the tetrahedral silicon lattice. The 



TABLE I: Character table for the double group Td- 
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(one-hole) ground-state is a ls-like fourfold degenerate 
state that belongs to the Tg representation of the tetra- 
hedral double group Td (for the nomenclature of repre- 
sentations used, see Table HJ. The bound hole has total 
angular momentum j = | . The single-hole wave function 
is the product of a Is hydrogen-like envelope function and 
a band-like function. Due to spin-orbit interaction in sil- 
icon, the j = \ valence band is split off by ~ 43 meV 
[l| and does not need to be considered in first order. A 
magnetic field completely lifts the fourfold degeneracy 
and the lowest order Zeeman effect of the Tg state is lin- 
ear. 

As far as symmetry is concerned, the B + -state is sim- 
ilar to neutral group-II acceptors in a tetrahedral lat- 
tice, which are well-studied (e.g. Ref. and references 
therein). Coupling two j = I (r 8 ) holes bound to a sin- 
gle nucleus gives rise to a six-fold degenerate state, be- 
cause due to the Pauli-principle only the tmta-symmetric 
part of Tg x Tg must be taken into account. This can be 
reduced to its components as {rg x Tg} = Ti + T3 + T5. 
Interaction between the two holes can split the state into 
a non-degenerate T\ state with total angular momentum 
J = and a fivefold T 3 + T 5 state carrying J = 2. 

Detailed quantitative calculations, which arc necessary 
to establish the ordering and splitting of these levels, are 
very difficult to carry out, because of the many compli- 
cated physical effects that must be taken into account 
(valence band structure, crystal field, Jahn- Teller-effect, 
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etc.). Hund's rule, well-known from atomic physics, pre- 
dicts that the more symmetric r 4 state has a higher en- 
ergy than the T 3 + T 5 state, such that the latter is the 
ground state. The same conclusion was drawn from nu- 
merical calculations based on effective mass theory 
However, it has been shown that a dynamic Jahn- Teller 
effect can provide a mechanism to reverse the ordering of 
the levels 111] , leading to a r 4 ground state. This has 
in fact been observed in several neutral double acceptors. 

Very little experimental work on B + has been done. 
The binding energy of the second hole in an isolated B + - 
state has been measured in phonon-induced conductivity 
(PIC) measurements and photoconductivity experi- 
ments ^3 . It is slightly below 2 meV. Stress-dependence 
has been investigated with the same techniques 0, El 
and in one case the results were explained as evidence 
for a stress-induced ground state splitting [l^. However, 
interpretation of the conductivity data is non-trivial, be- 
cause only levels which are very close to either the ground 
state of B + or the valence band edge can be observed 
with these techniques. Similar experiments in a mag- 
netic field 0| showed a linear increase of the binding 
energy, which was ascribed to Landau level formation in 
the valence band. In these experiments, no additional 
shift or splitting was resolved. 

Recent transport experiments in Si resonant tunneling 
devices provide a way to directly observe the magnetic 
field dependence of the B + state Q. These experiments 
showed a super linear shift of the ground state towards 
the valence band (Fig. |2J) . Neither a ground state level- 
splitting nor bound excited states were observed. 

DOUBLE ACCEPTORS IN A MAGNETIC FIELD 

Here, we present a group theoretical study to quali- 
tatively analyze the magnetic field behavior of isolated 
acceptors binding two holes in a tetrahedral semiconduc- 
tor for various possible states. This analysis is not only 
applicable to neutral group-II acceptors (e.g. Zn in Ge), 
but also to group-Ill acceptors binding an extra hole and 
singly ionized group-I acceptors (e.g. Cu~ in Ge). After 
this general part, we return to the specific situation of 
B+. 

We subsequently consider various possible two-hole 
levels and analyze their behavior in a magnetic field using 
perturbation theory. All such levels transform accord- 
ing to single- valued representations of X^, as shown in 
the overview in Table |n| We assume that the Coulomb 
force and spin-orbit interaction between the holes is suf- 
ficiently strong to split the levels into their irreducible 
components. Because of its possible importance for B + , 
we also consider the T 3 + T 5 level. In all cases it is as- 
sumed that the level under consideration is well separated 
from neighboring levels. 

Furthermore, we briefly address the analogue of the 



TABLE II: Overview of possible two-hole states arising from 
products of two single hole states and their reduction to ir- 
reducible representations of Td- For states originating from 
two equivalent single hole states (first two lines in the table), 
the Pauli principle allows only the antisymmetric part to be 
considered. 
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TABLE III: Character table for the double group S4 (lo = 

e ln/4 ). 
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central field approximation in atomic physics, where it 
is assumed that each of the two holes moves in the field 
of the negative ionized acceptor core and the averaged 
effective potential due to the other hole. In this approx- 
imation, the symmetry of the field in which each hole 
moves is unaffected by the presence of the second hole. 
This method is known to give a good description for some 
group-II acceptors in Si and Ge 0. 

When a magnetic field 13 is applied, new terms are 
introduced in the Hamiltonian of the holes, as given by 



TABLE IV: Character table for the double groups C3 (left; 
lo = e lw/3 ) and C lh (right). 
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TABLE V: Reduction of the representations of T d when a 
magnetic field is applied along a (100), (111), or (110) di- 
rection of the tetrahedral lattice. From this table it can be 
deduced how the double acceptor levels split in a magnetic 
field. 
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known phenomenological parameters, in terms of which 
the level splitting can be expressed. This approach is 
especially advantageous in the present situation, where 
both the values of the parameters occurring in Hz and 
the unperturbed wave functions are not (exactly) known. 



LINEAR ZEEMAN EFFECT 

In this section, we will investigate the first order Zee- 
man effect of all the double acceptor levels mentioned 
before. 



the Zeeman-Hamiltonian 
Hz = -ML + 2S)-© 

-lrn^ 2 B {(rl+rl)B 2 -[(r 1+ r 2 ).B} 2 }, 

where m* is the hole effective mass and L and S are the 
total orbital and spin angular momenta in units of H. The 
quantity L + 2S is the total static magnetic moment of 
the system. Moreover, \i& is the Bohr magneton and 
is the position vector of the i-th hole. 

The symmetry group of the Zeeman Hamiltonian TLz 
is Caoh- Unless B is directed along one of the main 
crystallographic axes, the symmetry group of the total 
Hamiltonian H = Ho + Hz reduces to the trivial group. 
When B is parallel to a (100), (111) or (110) direction 
in the crystal, the symmetry group of the total Hamilto- 
nian reduces from Td to S4, C3 or C\h, respectively. The 
relevant character tables are given in Table lllll and II VI 
Because all resulting groups are Abelian (commutative), 
it follows that the application of a magnetic field com- 
pletely removes the degeneracy of all levels 0] . The way 
in which the I\ levels exactly split in a magnetic field is 
presented in Table [V] 

To deduce the magnetic field induced splitting of the 
levels, we employ first order degenerate perturbation the- 
ory. As mentioned before, it is assumed that the sep- 
aration of the levels is large compared to the splitting- 
caused by the field, so only the subspace of Hilbert space 
connected to the level under consideration needs to be 
taken into account. Given a set of basis functions \i) for 
a particular level, we find the corresponding sub-matrix 
(i\Hz\j) of Tiz and diagonalize it to obtain the splitting 
as a function of B. 

Instead of trying to calculate matrix elements from TLz 
(after choosing a suitable set basis functions) it is much 
more convenient to use the well-established approach of 
constructing an effective Zeeman Hamiltonian 15]. This 
comprises the construction of a matrix of the required 
size, exploiting required symmetries to find vanishing el- 
ements and relations between elements. The result is a 
matrix that usually depends on a small number of un- 



The Tj-levels 

Because r 4 occurs in neither of the anti-symmetric di- 
rect products {Ti x Ti} (i = 1 ... 3), the effective Hamil- 
tonian matrix 7i c ff,im vanishes identically for the three 
levels IV Hence, none of these levels experiences a linear 
Zeeman effect. 

The linear part of the effective Zeeman Hamiltonian 
for a T4 or T 5 level is given by 0, Hif 

H e ff,lin = t>>Bg{B x J x + ByJy + B Z J Z ). 

Here, hb is the Bohr-magneton and J x , J y and J z are 
matrix representations of the components of the angu- 
lar momentum operator with respect to some conve- 
nient basis. The components J a (a = x, y, z) trans- 
form according to the r 4 representation of Td- Because 
{T 4 x T 4 } = {r 5 x r 5 } = r 4 , the r 4 and r 5 level do 
have a linear Zeeman effect. Calculating the eigenvalues 
of the matrix [W e ff,iin]i (* = 4,5) yields 



-HbgB. 

The eigenvalues are independent of the direction of the 
magnetic field and hence give rise to an isotropic split- 
ting. 



The T 3 + T 5 level 

The situation where the zero-field splitting of the 
T3 + T$ level is small compared to the Zeeman energy 
must be dealt with separately. Because T 3 x T 5 = r 4 + r 5 
contains T 4 , there are non vanishing cross-terms in the 
linear effective Zeeman Hamiltonian for a I^Ts-level. 
Therefore, such a level will have a linear Zeeman shift dif- 
ferent from that of the individual T3 and T$ levels. The 
Hamiltonian sub-matrix for the T3 + I^ level is given by 
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V a 35 B z J 

where 035 is assumed to be real. This 5 x 5-matrix is 
given with respect to a basis with two T3 wave functions 
and three Ts-wave functions. Only the matrix-elements 
connecting T 3 functions to T5 functions are shown ex- 
plicitly. The upper left and bottom right parts are the 
matrices of the individual T 3 and T 5 level, respectively, 
as given in the previous subsection. 

From this matrix, we determine the eigenvalues for 3 
parallel to the main crystallographic directions. For 3 \\ 
(100), so B x = B, B y = B z = 0, we find 

' 

AE = { ±^ B9 B 
±a 35 B. 

For 3 || (111), so B x =B y =B z = B/V%, we have 



where a\ is a phenomenological parameter. The simple 
conclusion is that a Ti level will experience a quadratic 
shift, independent of the direction of the magnetic field: 
AE = a\B 2 . From this purely symmetry-based analysis, 
( conclusions can be drawn neither about the magnitude 
of a\ nor about its sign (that is, whether the state is 
diamagnetic or paramagnetic). Because r2 X T2 = Ti, a 
similar expression holds for a T2 level. 



The T 3 level 

For a T3 level, the effective Hamiltonian contains two 
unknown parameters and is given by 0, Il7| 



AE 



Finally for 3 \\ 
it is found that 



±^ B gB±^fi 2 B g 2 + 2a 2 5 . 
(110), so B x = By = B/V2 and B z = 0, 



AE = 





±U 35 B 



±\B^2^ 



3a§ 5 - 



We conclude that there is indeed a linear Zeeman effect 
in the T 3 + T 5 level and the size of the effect is dependent 
of the direction of the field with respect to the crystal. 



n, 



err ,quad 



a 3 B 2 + b 3 [-(2B 2 -B 2 x -B 2 )a x 



V3(B 2 x ~B 2 y )a y ], 



where a x and a y are Pauli spin matrices and a 3 and 63 
are phenomenological parameters. When 3 || (100) the 
eigenvalues are given by 

AE = (a 3 ± 2b 3 )B 2 . 

This is a symmetric quadratic splitting superimposed on 
a quadratic shift. When 3 || (111), there is only one 
eigenvalue 

AE = a 3 B 2 , 

meaning that there is no splitting in second order and 
the quadratic shift is the same as for 3 || (100). Finally, 
for 3 || (110), we find the eigenvalues 

AE = (a 3 ± b 3 )B 2 . 

The Zeeman effect for this field direction is similar to 
3 || (100), but the splitting is twice as small. 



QUADRATIC ZEEMAN EFFECT 

For some of the levels we will also give a second or- 
der approach, using the quadratic part of the effective 
Hamiltonian 7Y c ff !qua d. Note that H e ff,quad contains both 
a second order approach to the linear part of the original 
TLz an d a first order approach to the quadratic part of 
the original TLz- 

The Ti and T 2 levels 

For the Ti level, the effective quadratic Zeeman Hamil- 
tonian contains only one term and is straightforwardly 
given by 

W e ff,q U ad = a,iB 2 , 



The T 4 and T 5 levels 

Because the symmetrized squares of and Ts satisfy 
[r 4 x T 4 ] = [r 5 x r 5 ], the results for the T 4 and T 5 lev- 
els are similar. For these two levels, the quadratic part 
of the effective Zeeman Hamiltonian has three unknown 
parameters a^, bi and Ci (i = 4, 5) and is given by 

^cff.quad = diB 2 

+h [3(B 2 X J 2 + B 2 y J 2 - B 2 J 2 ) - 2B 2 } 
+ Cl [B y B z {J y , J z } + B X B Z {J X1 J z } 

~]-B x By{J x: Jy\\ , 

where {A, B} = \{AB + BA) denotes the anti- 
commutator of A and B. We calculate the eigenval- 
ues of the full quadratic Hamiltonian matrix [H e ff,iin]i + 
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FIG. 1: Schematic overview of level splitting in a magnetic 
field along the main crystallographic axes in several kinds of 
two-hole acceptor levels. The figure illustrates the qualitative 
aspects of the splitting. The values of the parameter have 
been chosen to emphasize these features. 



[W c ff,quad]i for the three main crystallographic directions. 
For B || (100) we have 

fi B gB + (oj + b t )B 2 
AE = I ( ai - 2b l )B 2 

-imgB + (a, + bi)B 2 



in a magnetic field. The two-hole wave functions are anti- 
symmetrized products of one-hole wave functions and the 
energy levels are obtained by examining the various ways 
to put the two holes in the one-hole levels. 

We will present results for the case where both holes 
are put in a Tg level and B || (100) only. Similar results 
for the other types of levels and other directions of the 
field are easily obtained in an analogous way. 

For a magnetic field B || (100), the single hole Tg 
ground state is split into r 5 , T 6 , and r 8 -levels of 
£4 Q. Because holes are fermions, each of these non- 
degenerate levels can be occupied by at most one hole. 
By putting each of the two holes in a different level, this 
gives rise to six two- hole levels x Tg = T\ , r5 x Tj = T 3 , 

r 5 x r 8 = r 2 , r 6 x r 7 = r 2 , r 6 x r 8 =_r 4 , and 

IV x Tg = Ti, where all representations are of S4. 

The energy shifts of the single-hole levels have been 
determined experimentally Q. The shifts of the two- 
hole levels can be calculated as the sum of the shifts 
of the individual single hole levels from which they are 
composed. This results in a linear shift for each two-hole 
level, given by n B gB, with g = §g 3/2 + \g\/2 for T 2 , g = 
§33/2 - 5.91/2 for T 4 , g = for 2T U g = -|g 3 / 2 + 5P1/2 
for T 3 and g = —§.93/2 — 591/2 f° r ^2- The parameters 
<? 3 /2 and gi/2 are the g-factors for the single hole j = I 
and j = i-levels respectively. Experimental values for B° 
in Si are g 3 / 2 = 1.12 and g 1 / 2 = 1.04 0. In the above, a 
small overall shift is neglected. 



For B || (111) we have 



AE = I 



HsgB + (a, + \ci)B 2 



^ -[i B gB+ (ai + \ci)B 2 . 
And for B \\ (110) we find 



B^W + + \h)B 2 
+ (a 4 + \bt + \c t )B 2 
AE { (ai- \h - \ Cl )B 2 

-ByJ» 2 B g* + [\ Cl + \btfB 2 
+ (a 4 + \h + \ci)B 2 

It follows that in second order the spitting is no longer 
symmetric and isotropic for these levels. 

The Zeeman effect of the levels treated so far is 
schematically illustrated in Fig. ^ 



CENTRAL FIELD APPROXIMATION 

Finally we discuss the Zeeman effect for two-hole states 
in the central field approximation. In this approximation, 
we must start from the one-hole levels and their behavior 



APPLICATION TO B+ 

Only states arising from {Tg x Tg} (see Table [HJ are 
candidates for the B + ground state. These are T\ : T 3 , r 5 , 
T 3 + T 5 and the unsplit (central field) {T 8 x T 8 }. Each 
of these five possibilities for the B + ground state will 
be compared to existing experimental data. From the 
previous section, we conclude that all possible ground 
state levels behave qualitatively differently in a magnetic 
field. Therefore, it is in principle possible to determine 
the nature of the actual ground state of B + from the 
analysis of a sufficiently detailed experiment. Though 
this approach is hampered by the fact that the value 
of the parameters are not known, it is possible to draw 
conclusions based on the qualitative characteristics, such 
as linear or quadratic splitting/shift and the asymmetry 
of the splitting. 

We refer to our recent experiments reported in Ref. 
and summarize the main observations. The ground state 
energy shifts upwards (that is, in the direction of the 
valence band) and is therefore diamagnctic. The shift has 
both a linear and a quadratic component. The total shift 
amounts to 1 meV at a magnetic field of 14 T and was 
equal for the (100) and (110) directions (see Fig.|2J|. The 
width of the observed peak (full width at half maximum) 
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increased from 1.2 meV to 1.5 meV in the same magnetic 
field range. Within the experimental error (~ 0.2 meV), 
no splitting of the peak was detected |2l|. 

The experimentally observed super linear overall shift, 
independent from the direction of 3, best matches the be- 
havior of a Ti state, although this leaves the strong linear 
component in the measured magnetic field dependence 
unexplained. Therefore, we believe that the ground state 
of B + is indeed a Ti state. This hypothesis does imply 
that the observed linear component in the peak shift and 
the peak broadening are is due to other processes (e.g. 
the Stark effect) as already suggested in Ref. 0. 

The broadening in the observed peak is linear in the 
magnetic field and independent of its direction. There- 
fore, it cannot be explained as unresolved splitting of 
a T 3 level. A T 5 or T 3 + T 5 ground state would give 
rise to linear splitting (broadening), but no overall shift 
would be expected in first order. Moreover, the mag- 
nitude of the splitting in a r 3 + T§ level would depend 
on the magnetic field direction. Therefore, these possi- 
bilities are not consistent with the experimental obser- 
vations. Only when the parameter a 3 (05) would much 
larger than all other relevant parameters (that is a 3 <C 63 
or 05 <C &5,C5, Hbg/B), the magnetic field dependence of 
the T 3 (r 5 ) state would be similar to that of the Ti state. 
In that case, T 3 , T 5 and r 3 +T 5 states cannot be rejected 
as potential ground state symmetries for B + . 

The central field approach is unlikely to yield good re- 
sults for B + , for which the wave functions of the two holes 
are expected to overlap considerably (due to the small nu- 
clear charge). The peak splitting (or broadening, due to 
unresolved splitting) expected in this approach between 
the two T 2 levels would be given by 2^B(|g 3 / 2 + |5i/2)-B- 
Assuming the B°-values of the g-factors are valid here, 
this would amount to 3.6 meV for B = 14 T. This is 
much larger than the observed 0.3 meV increase of the 
FWHM of the measured resonance peak. Moreover, the 
1 meV shift observed in the experiment is much larger 
than the expected overall peak shift in this approach. 
Therefore, the description of the B + ground state in the 
central field approximation is not consistent with the ex- 
perimental observations. 

In summary, magnetic field dependent measurements 
indicate that the B + ground state is a non-degenerate Ti 
state. It would be interesting to have higher resolution 
spectroscopy data available, in order to exclude that the 
observed peak broadening is due to unresolved splitting. 
Knowledge of the B + ground state wave function would 
allow for obtaining quantitative information about the 
phenomenological parameters, which would be advanta- 
geous in the interpretation of experimental data. 
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FIG. 2: Magnetic field induced shift of the B + ground state 
deduced from electrical transport measurements Q. The an- 
gle between the magnetic field and the current is denoted by 9. 
The expression AE oc 13(1 + aB) was fit to the data, yielding 
a = 0.046 T _1 . 



CONCLUSIONS 

In conclusion, we have presented a general group the- 
oretical study of the magnetic field dependence of two- 
hole states of acceptors in tetrahedral semiconductors. 
We have used our results to analyze recent experimental 
observations. This analysis indicates that the B + ground 
state is most compatible with a Ti state. 
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